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Abstract 

We report on very strong evidence of the occurrence of power terms 
in OjjoM (p)' QCD running coupling constant in the MOM scheme, 
by analyzing non-perturbative measurements from the lattice three- 
gluon vertex between 2.0 and 10.0 GeV at zero flavor. While putting 
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forward the caveat that this definition of the coupling is a gauge de- 
pendent one, the general relevance of such an occurrence is discussed. 
We fit A^"'''* = 237 lb 3 '^-^ MeV in perfect agreement with the result 
obtained by the ALPHA group with a totally different method. The 
power correction to a^j[^{p) is fitted to (0.63 it 0.03 1 q'^^) GeV^/p^. 

1 Introduction 

The non-perturbative computation of the running QCD couphng constant 
as{p) follows a two-sided goal: the large energy matching to perturbative 
asymptotic QCD formula turns out to be a most direct method to predict Ay[§ 
from QCD first principles [jl| . On the other hand, the moderate or low energy 
behavior of as{p) is extremely instructive about non-perturbative properties 
of QCD. In this paper we restrict ourselves to high and intermediate energies 
and consider power corrections (~ l/p^) to the leading asymptotic behavior. 
As we shall see it turns out that there is no sharp separation between the 
asymptotic domain and the intermediate one. The power correction, beyond 
the lessons it contains by itself, greatly improves our asymptotic study and is 
never negligible up to ~ 10 GeV ! This surprising fact could only be revealed 
thanks to the high accuracy achieved in the present work. 

The asymptotic approach has been recently followed in Ref. p|. It is 
worth remarking that this matching procedure has also been developed for 
the two-point Green function to the same goal P], both matching programs 
leading to a consistent estimate of A^. The almost two-sigma discrepancy 
between the last estimate and the one obtained from Schrodinger functional 
methods ||^ seems to imply that some source of systematic uncertainty re- 
mains uncontrolled. Furthermore the careful study of the asymptotic be- 
havior carried out for the gluon propagator in Refs. stresses a danger: 
it is sometimes difficult to distinguish between the real asymptotic scaling 
and a behavior mimicking asymptoticity but with a certain effective "re- 
scaled" A parameter which differs significantly from the real asymptotic one. 
Consequently, and in spite of the agreement between both above-mentioned 
estimates of (matching two or three-point Green functions to perturba- 
tive formulae), we must inquire whether both are not biased by some sizable 
non-perturbative effect. 
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The operator product expansion (OPE) gives a standard procedure to 
parametrize non-perturbative QCD effects in terms of power corrections to 
perturbative results. In tliis framework, tlie powers involved in tlie expansion 
are expected to be uniquely fixed by the symmetries and the dimensions of the 
operators appearing in the product expansion. It should be noted that, due 
to the asymptotic nature of QCD perturbative expansions, power corrections 
are reshuffled between operators and coefficient functions in the OPE 

Since we work in Landau gauge, the gauge dependent local operator 
A^A^ is allowed in the OPE implying a dominant l/p^ power correction. 
Indeed sizable corrections are present as we shall show at length in the 
next section. 

In a less straightforward way our result is related to another hot issue 
(this is the spirit of the preliminary study in 0): the possible presence of 
terms in gauge m?7anant quantities such as Wilson loops P|. Since no gluon 
local gauge invariant operator of dimension less than 4 exists it is expected 
from OPE that the dominant power correction should be oc This 
picture has however recently been challenged [P-p!0[|. It was pointed out that 
power corrections which are not a priori expected from the OPE may in fact 
appear in the expansion of physical observables. Such terms may arise from 
(UV-subleading) power corrections to as{p), corresponding to non-analytical 
contributions to the /3-function. It is worth stressing the following. One 
knows that the perturbative analysis does not encode all the information 
on the coupling. Among all that is missed, a peculiar contribution to the 
coupling could result in a peculiar correction to physical observables. As a 
matter of fact, some evidence for an unexpected ^ contribution to the gluon 



condensate was obtained through lattice calculations in Ref. [g] (see also [|1T 
for an early evidence of such a contribution, although the perturbative series 
involved was not managed up to high orders). 

Let us insist: there is no direct relation between the corrections 
found in the present work to as{p) in a gauge dependent scheme and the 
power corrections advocated in the preceding paragraph to justify l/p^ cor- 
rections to gauge independent quantities. Nobody knows how to relate a 
gauge dependent scheme to a gauge independent one beyond the perturba- 
tive regime. Still, it may be that these two phenomena are not completely 
unrelated and the large l/p^ corrections found in the present work might 
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trigger some further thoughts along this hue. 

The paper is organized as follows: in Section ^ we explain the meaning 
of the lattice data, our strategy for the analysis and report the results. In 
Section R| we briefly review some theoretical arguments in support of power 

— A 2 

corrections to as(p), illustrating the special role of the p- term. Finally, in 
Section ^ we draw our conclusions. 

2 Lattice calculation of as from Green func- 
tions and power corrections 

2.1 as on the lattice 

Several methods for computing as{p) non-perturbatively on the lattice have 
been proposed in recent years |T^^T6|. In most cases, the goal of such studies 
is to obtain an accurate prediction for as{Mz), i.e. the running coupling at 
the Z peak, which is a fundamental parameter in the standard model. For 
this reason, lattice parameters are usually tuned so as to allow the computa- 
tion of as{p) at momentum scales of at least a few GeVs, where the two- loop 
asymptotic behavior is expected to dominate and power contributions are 
suppressed. 

As we shall see, however, non-perturbative power corrections cannot be 
neglected at energy scales as large as 10 GeV which is a sufficient reason 
to consider them in the fit. As a bonus the knowledge of these power cor- 
rections provides us with a physically significant quantity as argued in the 
introduction. 

For this purpose, the best method is one where one can measure as{p) in a 
wide range of momenta from a single Monte Carlo data set. Indeed, a narrow 
energy window does not allow to disentangle in a clear cut manner the power 
corrections from unknown higher perturbative orders, and these corrections 
can be mimicked by an effective Aqcd different from the asymptotic one. 

One method which fulfills the above criterion is the determination of the 
coupling from the renormalized lattice three-gluon vertex function [0,0,0]. 
This is achieved by evaluating two- and three-point off-shell Green's functions 
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of the gluon field in the Landau gauge, and imposing non-perturbative renor- 
mahzation conditions on them, for different values of the external momenta. 
By varying the renormalization scale p, one can determine as{p) for different 
momenta from a single simulation. Obviously the renormalization scale must 
be chosen in a range of lattice momenta such that both finite volume effects 
and discretization errors are under control. Both these constraints impose 
too strict limits on the energy range if only one lattice run is used. This is 
why the procedure used in Ref. combining several lattice simulations at 
different /3's is a necessity to get a larger range of lattice momenta. The use 
of different volumes will also help to control finite volume artifacts. 

The definition of the coupling corresponds to a momentum-subtraction 
renormalization scheme in continuum QCD |T^. It should be noted that 
in this scheme the coupling is a gauge- dependent quantity. As we already 
mentioned, one consequence of this fact is that corrections should be 
expected, based on OPE considerations. For full details of the method and 
the lattice calibration we refer the reader to Ref. 

2.2 Models for power corrections and construction of 
the data set 

In the present work we shall not address the general problem of defining an 
optimal analytic form for the coupling at all scales to which we could fit 
our data. For the purpose of our investigation, we shall compare the non- 
perturbative data for as with simple models obtained by adding a power 
correction term of the form l/p"^ to the perturbative formula at a given or- 
der. This amounts to a quite raw separation between a perturbative versus 
a non-perturbative contribution, the major problem of course being the pos- 
sible interplay between a description in terms of (non-perturbative) power 
corrections and our ignorance about higher orders of perturbation theory. 
As crude as it may be, our recipe allows for addressing this problem, as we 
shall see. 

In order to identify a momentum interval where our ansatz could fit the 
data, one should keep in mind that the momentum range should start well 
above the location of the perturbative Landau pole, but it should nonetheless 
include low scales where power corrections may be large up to large enough 
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momentum scales in order to be confident that the asymptotic regime {i.e. 
perturbative running) has become dominant. Our choice will be| 2.0 GeV 
< p < 9.6 GeV. It will turn out that in the whole range both the perturbative 
three-loop contribution and the non perturbative l/p"^ one contribute by a 
sizable amount, although obviously the former becomes dominant at larger 
scales. 

A data set which fulfills the above-mentioned requirements can be con- 
structed by aggregating data computed at different P values (/3 =6.0, 6.2, 
6.4, 6.8) on a 24^ lattice. The fact that such a data set can be assembled is 
a very good a posteriori check that the expected scaling in the lattice cutoff 
a takes place. 

On the other hand, the physical volumes for these simulations are very 
different from each other. The appropriate matching of the data proves that 
finite-size effects remain controlled. The pattern for these volume effects is 
clear from fig. |l|(a), where the whole set of data is plotted, including the 
points too much affected by finite volume artifacts and eventually rejected 
in our fits. As can be seen, the effects are negligible for large enough Lp, 
where L is the physical lattice length. Evidences for the same trend arise 
from the comparison of data obtained on two different lattices (16^, 24^) 
aX (3 = 6.8, shown in fig. |l](b). In practice we will take as infrared cut-off 
at each (3 the value Pmin such that, including points below this value, the 
Xdof increases dramatically. This criterion leads to: Pmm(6.0) = 2.0 GeV, 
Pmin{G-2) = 2.5 GeV, p™„(6.4) = 4.0 GeV and p™„(6.8) = 6.0 GeV. It 
corresponds to Lp > 24 ^ Reversely, when we vary the infrared cut-offs 
above the values just mentioned, the Xd.o.f stays stable, in the range 1.4 to 
1.6. Even more striking, the resulting value for A^ is very stable, varying 
no more than 2 MeV. Consequently the data set obtained with these cut-offs 
should be considered as IR safe and will be used in the following fits. 

At each j3 value one should of course also worry of data in the range 
of the larger values of momentum, which are affected by lattice artifacts of 
0{a?p'^) (UV discretization effects). This was taken care of by the sinus- 

■^The range's upper limit is determined by the condition a p < 7r/2, where a is the 
lattice spacing, applied to our 24^ lattice at f3 — 6.8 . 

''Almost all the data from our 16'' lattice at f3 — 6.8 turn out to be excluded by such a 
requirement, Lp > 24 leading for example to an infra-red cut-off equal to 9. GeV. This is 
why we will not use at all the volume 16*. 
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Figure 1: Evaluations of as from 24^ lattices at /? = 6.0,6.2,6.4,6.8 are respec- 
tively shown in plot (a) with black circles, white diamonds, black squares and 
white circles. In plot (b) black (white) circles correspond to as evaluations from 
a 24"^ (16"^) lattice at /3 = 6.8 . 



7 



improvement program which has aheady been described in [0, which ba- 
sically amounts to the substitution of the lattice momenta with their 
O(a^p^) analogues |sin(^). As already noticed in 0, a good rationale for 
this in our gauge-fixed situation is that the gauge fixing algorithm leads to 
the relation ^ sin(^) A^{p) = 0, while p^A^ does not vanishQ. It should be 
noticed that without this prescription the quality of almost any fit degener- 
ates. The relevant configurations were generated on a QHl Quadrics system 
based in Orsay (see p). 



2.3 Fitting data to our ansatz 

We now proceed to fitting the data set to our ansatz, which, we recall, is the 
addition of a term proportional to ^ to a given order of perturbation theory: 



(1) 



Working at three loop level, the perturbative expression for the running cou- 
pling constant, a^'^^^{p'^) requires to inverse either the unexpanded formula 
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or the expanded one 



A.AM(S )(l.g|=^5) (3) 

^One should keep in mind that we are imposing Landau gauge. 
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where A^'^-' denotes the conventional two loops formula: 



and A = 2pop2-4:pf > (for our MOM scheme). In the previous formula, the 
use of A, 5 and /?'s stands for the A parameter, the running coupling constant 
and beta function coefficients in the particular MOM renormalization scheme. 
From now on we will systematically convert A into Aj^ using 

Aj^ = exp(-70/66) A ~ 0.346 A (5) 

In Eqs. (0-^ the p"^ dependence of 5 has been omitted for simplicity. 

No analytical expression can exactly inverse neither unexpanded eq. (|^) 
nor expanded eq. @. The following formula gives an approximated solution 
to the inversion of eq. (|^): 



?^r„2^ _ i!L _ SttA log(t) 

Pot f3o i/3oty 

where t = log(p^/A^). On the other hand, an exact numerical inversion of 
Eq. (H) can be easily obtained and used to fit our data. Both ansatze, eq. (|^) 
and the numerical inversion of eq. (0), should only differ by perturbative 
contributions higher than three loops. Thus, we will fit with the two ansatze, 
the difference between these two predictions being considered as an estimate 
of the systematic uncertainty coming from the neglected perturbative orders. 
To make more direct the comparison with the Schrodinger functional estimate 
of A^ the central value for our prediction of Aj^ should be taken from 
the fit with the exact inversion of Eq. (^Q. This yields: 

®The ALPHA collaboration estimate: A^jg — 238(19) MeV comes from evaluating 
eq. (pf) for a value of as obtained from the lattice at very high momenta. 
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= 237 ± 3 MeV, c = 0.63 ± 0.03 GeV^ xVd-o.f. = 1.6 ; (7) 



in perfect agreement with the determination 0] which uses a totally different 
technique. Using eq. (|) we obtain Aj^ = 227(5) MeV and c = 0.50(6) GeVl 
Comparing the latter results with eq. (^ provides an estimate of the higher 
loop uncertainty of about 10 MeV for and 0.1 GeV^ for c. The size of 
the power correction will be discussed in sections 3.1 and 3.2. 



' ' ' I 




1.1 GeV 

Figure 2: Formula (||), with the values of and c given by eq. (|^, fits impres- 
sively the data set built as explained in the text 

Fig. 3 illustrates the effect of the power correction in a striking manner. 
The upper set of points shows A, converted to A^, computed through eq. (|^) 
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from as provided by the lattice at every value p. Scaling would imply a 
constancy of which is far from true. The lower set of points corresponds 
to the same formula applied to ^^(p^)^^**"^'^ — 0.63/p^, i.e. to what should 
fit the perturbative formula. Now the constancy of is very good within 
the statistical errors. It is now clear why the upper points show a trend to 
decrease: as the energy scale increases, the effect of the power corrections 
must decrease, and the upper points converge slowly towards the lower ones. 
The surprise is that above 9.0 GeV the merging of the two sets of points 
has not yet taken place contrarily to the general expectation that power 
corrections are negligible at such a scale. We will elaborate on this in the 
next section. 

We note that using eq. (^) one can draw the same conclusion: imposing 
c = (i.e. fitting to a pure three loop formula) there is no good fit on the 
whole range of momenta and the best that one can obtain on a restricted 
interval yields a value for Ay^ higher than expected (e.g. with respect to 0]) 
and a definitly worse Xd.o.f- 

We now proceed to address the question of the stability of our previous 
result with respect to the inclusion of the fourth loop. To fit the value of A, 
one would need /^s, the four loop coefficient in the perturbative expansion of 
the /3-function. This coefficient is unknown but one can pin down a reason- 
able range of variation. Let us consider the ratio &3/62 with 63 = jS^/lATc)^ 
and 62 = /52/(27r)^. This ratio is larger than one in the only scheme (MS) 
for which it is known. On the other hand we expect that, for a reasonably 
small value of the coupling, the perturbative expansion of the /3-function is 
still in the regime in which it seems to converge at four loops. As a working 
hypothesis, we then suppose that the contribution to the /3-function coming 
from jSs (that is, the one proportional to a^) is not too much larger than the 
one coming from /?2 at a typical value of a ~ 0.4, implying 63/62 < 2.5. Ac- 
tually we conservatively vary 63 from up to 562 in the following exploration 
of a large range of values for 63. 

First of all we try a pure four loop fit (that is without any power correc- 
tion). We observe that there is no good fit on the whole range of momenta. 
If one tries to add again a term proportional to ^ to the four loop pertur- 
bative expression, the following should be notecf Good fits are recovered 
either on the whole range of momenta for 63/62 ^ 1 or by discarding mo- 
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Figure 3: The estimates of computed from ag via eq. ^ are plotted as 
a function of the renormalization point momenta. Diamonds, black and white 
circles, and asterisks respectively correspond to evaluations from 24^ lattices at j3 = 
6.8,6.4,6.2,6.0. Upper points are computed from the values directly obtained 
from the lattice, while the lower points use af as defined in eq. with c given 
by eq. ^). 



menta below 3 GeV as 63/62 ^ 4. As far as the value of A is concerned, it is 
really stable when the fits are of good quality^} As for the coefficient of the 
non-perturbative term ^, it is less stable. Results are summarized in Tab. |1]. 
One can see how the "new player on the ground" , /^s, is strongly correlated to 
the coefficient of the power term. This coefficient is anyhow fully consistent 
with the value determined at three loop level, as long as 63/62 < 2 i.e. when 
the asymptotic behavior is not too dramatically perturbed by the four-loop 
contribution. 

^Even for 63/62 — 5, one could obtain a good quality fit (x^ = 1-8) with = 238, if 
momenta below 4 GeV are now discarded. 
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7.8 


299 


8.3 


294 


1.6 


237 


0.63 


1.6 


235 


0.67 


1 


7.1 


284 


6.0 


287 


1.8 


238 


0.56 


1.6 


235 


0.61 


2 


28 


259 


4.1 


279 


2.7 


229 


0.57 


1.6 


236 


0.53 


3 


104 


227 


3.0 


270 


4.7 


215 


0.65 


1.6 


238 


0.45 


4 


145 


215 


2.9 


261 


7.3 


202 


0.75 


1.8 


237 


0.37 


5 


157 


209 


4.4 


252 


10.2 


190 


0.84 


2.4 


231 


0.37 



Table 1: A collection of fitted parameters obtained by imposing different values for 
the ratio 63/62 defined in the text. "Whole" refers to the whole energy window (2 
GeV < fi <10 GeV) and "> 3" to a momentum range above 3 GeV. 



Again, much the same holds when fitting to a formula for function 
of momentum (like eq. (^) at the four loop level. In order to trust a four loop 
more than power corrections one would need both to discard lower momenta 
and to accept excedingly large values for 63. 

Of course, the perturbative knowledge of P3 coefficient is unavoidable to 
get total confidence on our results versus higher loop orders inclusion. One 
should take this last analysis only as a preliminary check. Still it is another 
indication that things are working pretty well (that is, consistently) with 
respect to our theoretical prejudice. In any case, the value obtained for Af^ 
is almost insensitive to P^. 

We collect all the hints from these many counterproofs. To estimate the 
systematic errors we use two methods: first we compare the exact inversion 
of eq. (^ and the use of eq. on our results, second we use the results on 
the whole energy range in table with a reasonable x^- We take the three 
loop result with formula (^ as our central value 

Aj^ = 237 ± 3 til MeV, c = 0.63 ± 0.03 t GeV^ (8) 

The present analysis has improved over the one presented in [|l],3 by 
using a larger data set which provides a wider momentum window and taking 
into account a power correction. If one tried to repeat the fits only in the 
range of the lower momenta as we did in []T], , there would be no clear cut 
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indication for power corrections and the effective value for A would turn 
out to be higher. What the upper momenta data really do is to single out 
the asymptotic value for A, while deviation from asymptotia in the lower 
momenta data asks for power correction rather than higher loops. 



3 Discussion on the ^ corrections to as{p) 

Power corrections to as{p) can be shown to arise naturally in many physical 



schemes |]18|,|19[. The occurrence of such corrections cannot be excluded a 
priori in any renormalization scheme. Even more so, as previously stated, 
in a gauge dependent renormalization scheme as the MOM discussed here. 
Clearly, the non-perturbative nature of such effects makes it very hard to 
assess their dependence on the renormalization scheme, which is only very 
weakly constrained by the general properties of the theory. 

As discussed in the following, several arguments have been put forward in 
the past to suggest that a most likely candidate for a leading power correction 
to as{p) would be the same term of order A^/p^ we found. Furthermore, it 
is worth emphasizing that this does not contradict the OPE expectation for 
a gauge dependent quantity. 



3.1 Static quark potential and confinement 

Consider the interaction of two heavy quarks in the static limit (for a more 
detailed discussion see [0). In the one-gluon-exchange approximation, the 
static coulombic potential V{r) can be written as 

V{r) ^ as I d'k^-^. (9) 

Using standard arguments of renormalon analysis, one may consider a 
generalization of (^ obtained by replacing a.^ with a running coupling: 

V{r) cx / d'kaMm^-^- (10) 
J \k\ 
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The presence in as{k'^) of a power correction term of the form c/k"^ would 
generate a hnear confining piece Kr in the potential V{r). It would of course 
be totally unjustified to take seriously our power correction in order to de- 
rive the linear potential. If we nevertheless perform this exercise to have a 
feeling of the scales, we get = 2/3c ~ 0.4 GeV^, while the usual string 
tension is around 0.2 GeV^. Note that a standard renormalon analysis of 



p!0D (see for the details) reveals contributions to the potential contain- 



ing various powers of r, but a linear contribution is missing. This is a typical 
result of renormalon analysis: renormalons can miss important pieces of non- 
perturbative information. 



3.2 An estimate from another lattice method 

The lattice community has been made aware for some time of the possi- 
bility of non-perturbative contributions to the running coupling; for a clear 
discussion see ||21| . Consider the "force" definition of the running coupling: 



"™-<« = r'^ W = ;). (11) 

where again V{r) represents the static interquark potential. By keeping into 
account the string tension contribution to V{r), which can be measured in 
lattice simulations, one obtains a l/Q^ contribution, whose order of mag- 
nitude is given by the string tension itself. Ironically, this term has been 
mainly considered as a sort of ambiguity, resulting in an indetermination in 
the value of a{Q) at a given scale. From a different point of view, such a 
term could be interpreted as a clue for the existence of a contribution, 
and it also provides an estimate for the expected order of magnitude of it, 
at least in one (physically sound) scheme. The same naive exercise than in 
the preceding subsection leads from eq. (|11]) to i^' = 4/3c ~ 0.8 GeV^. 

In order to make a deeper contact with what we are actually studying, 
one should of course keep in mind that all the preceding arguments would 
be physically sounder in the Coulomb gauge in which the notion of force 
has a clearer meaning. In Landau gauge such a picture is far from clear. 
Furthermore, it is known that the interquark potential becomes close to linear 
at distances larger than 0.5 fm, which corresponds to a momentum smaller 
than 0.4 GeV, while our l/p^ fit only apply above 2.0 GeV i.e. at distances 



15 



smaller than 0.1 fm. Still, a couple of considerations are in order at this point. 
We note that the rough order of magnitude of the power term we found is the 
same as inferred from the simple arguments from the static potential (some 
10~^ GeV^), even if we are in a different (UV) regime. While this could turn 
out to be accidental, it is nevertheless intriguing to think of some sort of 
relation. In any case, one should nevertheless note that the power correction 
that we have obtained is rather large with respect to the common wisdom 
of non-perturbative effects being negligible at scales such as 10.0 GeV. A 
further study of the relation between the power correction and the confining 
potential is clearly needed. 

3.3 On the relation with the lattice gluon condensate 



Having just stressed that the contribution of the power correction is rather 
large also at a scale such as 10.0 GeV, this is a good point to go back to 
the argument referring to the unexpected result of P] . As already mentioned 
in the introduction, non-perturbative contribution to the running coupling 
can be advocated in order to explain the A^/Q^ contribution to the gluon 
condensate. The argument runs as follows (see [Q). From general arguments 
one expects the condensate W to be written in the form 



which is based on the fact that this condensate has dimension four and is 
renormalization group invariant, so that the function /(p^/A^), which is in- 
dependent on Q (for large Q), can be expressed as a function of a running 
coupling. This leads to consider the contribution coming from the large fre- 
quency part of the integral 



puzzle. 




(12) 




(13) 
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in which the function f{p^/h?) is taken proportional to the running couphngQ. 
By taking into account the perturbative running couphng the IR renormalon 
contribution can be obtained (see for details). Again, one can consider 
also contributions coming from a non-perturbative correction to the coupling 
of the form c/p^. From the UV limit of integration one then obtains a 
A^/Q^ contribution to W . Let us insist: this is a contribution coming from 
a coupling with a c/p^ correction in the UV region. While stressing again 
that one can draw no definite conclusion from the following consideration, 
still it is worth noting that our scheme provides an example of a coupling in 
which a c/p^ correction is not negligible even at 10.0 GeV. For a discussion 



of possible scenarios for the result of [S] see also p2 



3.4 Landau pole and analyticity. 

It is well known that perturbative QCD formulae for the running of in- 
evitably contain singularities, which are often referred to as the Landau pole. 
The details of the analytical structure depend on the order at which the (3- 
function is truncated and on the particular solution chosen. The existence of 
an interplay between the analytical structure of the perturbative solution and 
the structure of non-perturbative effects has been advocated for a long time 
To illustrate this idea, consider the one-loop formula for the running 
coupling asip): 

aJp'^) = (14) 

bo log(|l) ^ ^ 

Here the singularity is a simple pole, which can be removed if one redefines 
as{p) according to the following prescription: 



1 A^ 

bo log{^) &o(A2-p2) 

where a power correction of the asymptotic form ^ appears. However, the 
sign of such a correction is the opposite of what one would expect from the 
results of and from the results in Section || (although the absolute value 

^Note that higher powers would be subleading both for renormalons and for power 
corrections. 
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is of the right order), so that in the end one could envisage a more general 
formula for the regularized coupling: 

. 1 A2 

bo log(^) bo{A'-p') >2 y J 

The message from (^) is that the coefficient of the power correction is not 
constrained by the mere cancellation of the pole. 

At higher perturbative orders one encounters multiple singularities, which 
include an unphysical cut. There are several ways to regularize them. In 
particular, the method discussed in combines a spectral-representation 
approach with the Renormalization Group. The method was originally for- 
mulated for QED, but it has recently been extended to the QCD case ||23 . 



4 Conclusions 

We have studied the strong coupling constant estimated non perturbatively 
on the lattice from Green functions in the Landau gauge using the MOM 
scheme. This has been performed with a large statistics of 1000 field config- 
urations per run and running at /3 = 6.0, 6.2, 6.4, 6.8. Finite volume effects 
as well as finite lattice spacing effects have been carefully controlled. 

We have parametrized the momentum dependence of as using the three- 
loop perturbative formula plus a c/p"^ term. We have obtain a good fit in the 
energy range from 2.0 GeV up to 10.0 GeV. As a result of this study we find 



A^js = 237 ± 3 tiQ MeV, c = 0.63 ± 0.03 t oa3 GeV^, (17) 

^Ms agrees perfectly well with the result of and the existence of sizable 
power corrections is convincingly established. The stability of our fit has been 
extensively checked. 

The power correction turns out to be rather large, providing a 3% cor- 
rection on as at 10.0 GeV, i.e. a 20 % correction on A. 

Having gathered evidences for a particular scheme, one needs to address 
the issue of the general relevance of such a finding in the spirit of the dis- 
cussion of section H, assessing the scheme dependence of our results. As 
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already discussed, the non-perturbative nature of power corrections makes it 
very hard to formulate any theoretical procedure to estimate the impact of 
scheme dependence. The best one can do at this stage is to consider different 
renormalization schemes and definitions of the coupling and gather numer- 
ical evidence and formal arguments supporting power corrections to as{p)- 
In this way, scheme-independent features may eventually be identified. For 
example, on the basis of our results, we note the following: 

• Theoretical arguments suggest l/p"^ corrections both for the coupling 
as defined from the static potential and for the one obtained from the 
three-gluon vertex. The arguments for the former case were outlined 
in Sections 3.1 and 3.2. As far as the coupling from the three-gluon 
vertex is concerned, 1 /p^ corrections appear in an OPE analysis if one 
keeps into account the fact that for such a gauge dependent coupling 
dimension 2 condensates are expected. 

• In the static potential case, the theoretical arguments also provide an 
estimate for the order of magnitude of the coefficient of the 1 /p^ correc- 
tion while in the three-gluon vertex case the OPE arguments do not, 
suggesting instead that it may depend on the gauge. Nevertheless , 
the order of magnitude of our numerical result in the Landau gauge is 
roughly the same as the one from the static potential case. This calls 
for further investigation, which may be performed for example by at- 
tempting a similar calculation in a different gauge and particularly in 
the Coulomb gauge in which the static potential is naturally defined. 

This issue of scheme dependence could be the focus of possible future 
work. 
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